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I.  Introduction 

High  average  power  lasers,  e.g.,  free  electron  lasers  (FELs)  and  solid-state  lasers 
(including  fiber  lasers)  are  prime  candidates  for  efficient,  directed  energy  applications. 
These  include  laser  power  beaming  and  laser  weapons,  requiring  multi-kWs  of  CW 
power  operating  in  the  IR  regime. 

High-intensity  lasers  operate  in  a  different  regime,  e.g.,  peak  powers  of  ~  1012  - 
1015  W,  pulse  lengths  of  -  10"12  -  10"14  sec,  intensities  of  ~  1014  -  1023  W/cm2,  and 
repetition  rates  ranging  from  10  -  10  Hz  with  average  powers  of  >  10’s  W.  These 
lasers  are  used  in  high-field  physics  research  and  have  numerous  potential  applications. 
The  physical  processes  associated  with  ultrashort  pulse  lasers  (USPL)  interactions 
include:  photo/collisional  ionization,  Kerr  and  rotational  Raman  effects,  self  phase 
modulation,  filamentation,  group  velocity  dispersion,  optical  shocks,  frequency  chirping, 
etc.  We  discuss  these  interrelated  physical  processes  and  some  of  the  applications  of 
these  lasers. 

One  of  the  important  topics  to  be  discussed  in  connection  with  both  high-power 
and  high-intensity  lasers  is  propagation  in  a  turbulent  atmosphere.  Laser  propagation  in 
atmospheric  turbulence  can  results  in  beam  centroid  wander,  spreading  and  intensity 
scintillation.  A  phase  conjugation  technique  to  mitigate  the  effects  of  atmospheric 
turbulence  is  described. 

In  Sec.  II  we  discuss  free  electron  lasers  and  related  topics  including  optical 
guiding,  efficiency  enhancement,  electron  beam  quality,  coherent  and  incoherent  x-ray 
FELs  and  backward  Raman  amplification.  In  Sec.  Ill  we  discuss  ultrashort  pulse  laser 
interactions  and  related  topics,  including  atmospheric  propagation  of  USPLs  and  laser 
driven  electron  acceleration  in  tapered  plasma  channels.  In  Sec. IV  we  discuss  the  role  of 
atmospheric  turbulence  on  the  propagation  of  laser  beams  and  also  briefly  discuss  optical 
phase  conjugation  as  a  means  to  mitigate  the  effects  of  turbulence.  In  Sec.  V  we  discuss 
various  applications  of  high-power,  high-intensity  lasers  including  remote  detection  of 
radioactive  material  using  electromagnetic  signatures,  atmospheric  lasing  of  N2 
molecules.  Concluding  remarks  are  presented  in  Sec.  VI. 


2 


II.  Free  Electron  Lasers 

Free  electron  lasers  can  generate  high-power  coherent  radiation  from  the 
microwave  to  the  x-ray  regime  [1-3].  The  FEL  mechanism  is  based  on  classical  physics 
in  which  an  electron  beam  induces  radiation  as  it  traverses  a  periodic  magnetic  field 
(wiggler)  [4-13].  The  FEL  can  operate  in  the  low  gain  (interference,  oscillator  regime) 

[4]  or  a  high  gain  (exponential,  amplifier  regime)  [5].  One  of  the  important 
characteristics  of  FELs  is  that  the  generated  radiation  beam  can  be  guided  by  the  electron 
beam,  making  possible  extended  interaction  lengths.  Also  FEL  efficiency  can  be 
enhanced  by  various  methods  such  as  spatially  tapering  the  wiggler  field.  Free  electron 
lasers  are  generally  large,  complex  and  expensive. 

A  schematic  of  an  FEL  is  shown  in  Fig.  1.  In  this  configuration,  electron  bunches 
from  a  radio  frequency  (RF)  linac  propagate  through  a  wiggler  and  radiation  is  amplified 
by  the  electron-wiggler  interaction.  The  spent  electron  beam  is  deflected  either  into  a 
beam  dump  or  re-circulated. 


relativistic 
electron  beam 
(RF  Linac) 


wiggler 


optical 

beam 


Figure  1 .  FEL  amplifier  configuration 


Figure  2  is  a  schematic  of  an  FEL  employing  a  planar  wiggler  with  period  that  is 
polarized  along  the  x  axis.  An  electron  beam  with  relativistic  factor  y  propagates  along 
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the  z-axis  and  wiggles  in  the  y-z  plane.  Radiation  with  wavelength  k  that  is  polarized  in 
the  x  direction,  propagates  in  the  z  direction  and  is  amplified  provided  A  «  /  2  y2  . 
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Figure  2.  An  electron  beam  interacts  with  a  wiggler  field  and  amplifies  radiation 
propagating  along  the  z-axis.  The  electron  wiggle  motion  is  in  the  y-z  plane. 


In  the  FEL  a  relativistic  electron  beam  propagates  through  a  combination  of  a 
wiggler  field  B„,  =  Bw  cos  (kw  z)ex ,  and  an  optical  radiation  field 

Es  =  Esin(kz  -  cot)  ey,  where  kw  —  2  71  /  is  the  wiggler  wavelength, 

k  »  col c  =  2x11,  X  is  the  radiation  wavelength  and  e  t ,  e  denote  transverse  unit 

vectors.  The  essential  feature  of  the  FEL  mechanism  is  that  the  beating  of  the  fields 
causes  the  electrons  to  undergo  axial  bunching  generating  a  beat  wave.  The  axial  force 
on  the  electrons  is  proportional  to 

(q!c)(\wx Br  +  vsxBw)z  ~  sin[(&  +  £w)z  -  £»f]ez, where  vw(vx)  isthe 

transverse  electron  velocity  due  to  the  wiggler  (radiation)  field.  This  beat  wave  is 
referred  to  as  the  ponderomotive  wave.  The  phase  velocity  of  the  ponderomotive  wave 
co/(k+kw )  <  c  is  slightly  less  than  the  speed  of  light  and  comparable  to  the  electron 

axial  velocity  v  _ .  The  relativistic  electrons  can  therefore  transfer  kinetic  energy  to  the 
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radiation  field,  resulting  in  amplification.  The  electrons  eventually  become  trapped  in  the 
ponderomotive  wave  and  the  radiation  field  saturates. 

Greater  insight  into  the  FEL  mechanism  is  provided  by  the  one  dimensional 
dispersion  relation.  The  electron  density  wave,  together  with  the  electron  wiggle  motion, 
generate  a  transverse  current  having  the  proper  phase,  i.e.,  spatial  -  temporal,  dependence 
to  amplify  the  radiation  field.  The  amplified  radiation  in  turn  enhances  the  electron 
density  wave  further  amplifying  the  radiation.  Considering  spatial  variations  in  the  z 
direction  only  the  FEL  dispersion  relation  is  [2] 

[k2  —  co2  / c2  )(o-vz(fc  +  kw))2  -co2p  / /3  )=  -4(o)2p  / r)a^  k2 ,  (1) 

where  aw  =  qBw  /me kw  is  the  normalized  wiggler  amplitude  and  cob  is  the  beam 

plasma  frequency.  The  dispersion  relation  describes  how  the  wiggler  field  couples  the 
electromagnetic  modes  (first  factor  on  the  left  hand  side)  to  the  electron  beam  modes 
(second  factor  on  the  left  hand  side). 

a)  Optical  Guiding  in  FELs 

In  the  FEL  the  wavefronts  of  the  amplified  radiation  are  such  that  the  radiation 
tends  to  focus  toward  the  electron  beam  [13].  Solving  the  dispersion  relation,  Eq.  (1),  the 
radiation  wavenumber  is  found  to  be  given  by  k  =  co/c  +  Ak  -  iT ,  where  the 
wavenumber  shift  M  >  0  when  the  growth  rate  is  T  >  0  (amplification).  Since  the 
refractive  index  associated  with  the  radiation  is  n  =  l  +  cAk/a>  >  1,  the  radiation 
undergoes  focusing  [13].  When  focusing  is  balanced  by  diffraction  the  radiation  beam  is 
optically  guided  though  the  interaction  region  as  shown  in  Fig.  3.  Optical  guiding  is 
important  since  the  interaction  length  in  the  wiggler  is  generally  long  and  diffraction 
would,  in  the  absence  of  guiding,  spread  the  radiation,  terminating  the  FEL  interaction. 

Figure  3.  Radiation  beam  is 
guided  (Ak  >  0)  by  electron 
beam  in  an  FEL  amplifier. 
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b)  Efficiency  Enhancement 

Amplification  of  the  radiation  is  at  the  expense  of  the  electron  beam  kinetic 
energy.  As  the  electrons  slow  down  the  PEL  interaction  coupling  is  reduced  and 
efficiency  decreases.  However,  by  spatially  tapering  the  wiggler  field  it  is  possible  to 
reduce  the  phase  velocity  of  the  trapping  (ponderomotive)  wave  to  maintain  synchronism 
with  the  electrons  [11].  This  is  one  of  a  number  of  methods  to  enhance  the  PEL 
efficiency.  Figure  4  is  a  plot  of  the  output  power  versus  interaction  length  obtained  from 
a  simulation  using  the  code  GENESIS  [14].  The  electron  and  wiggler  parameters  on 
shown  in  the  figure.  The  solid  curve  is  for  an  untapered  wiggler  resulting  in  ~  80  MW  of 
output  power  while  the  dashed  curve  corresponds  to  a  tapered  case  resulting  in  ~  500 
MW. 
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Figure  4.  GENESIS  simulation  of  FEL  power  versus  interaction  length  without  (solid 
curve)  and  with  (dashed  curve)  tapering.  The  insets  show  the  electron  phase  space 
(energy  vs.  ponderomotive  phase)  at  various  positions. 


Theoretically  and  experimentally  it  has  been  demonstrated  that  the  FEL  mechanism  can 
be  reversed  resulting  in  a  laser  driven  electron  accelerator  [15]. 
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c)  Electron  Beam  Quality 

In  the  FEL,  the  gain,  growth  rate  and  efficiency  are  sensitively  dependent  on  the 
axial  velocity  spread  A/?,  of  the  electrons.  The  electron  beam  quality  is  measured  in 
terms  of  emittance  and  energy  spread,  the  origins  of  which  include  i)  temperature  of  the 
cathode,  ii)  roughness  of  the  emitting  surface,  iii)  non-uniformity  of  the  emission  from 
the  cathode  surface,  iv)  nonlinearity  of  electric  and  magnetic  fields,  and  v)  self-fields. 

The  electron  beam  quality  requirement  for  the  FEL  is  that  A/7,  be  small  compared  to  the 
efficiency  A  /  2 Leff  at  the  entrance  to  the  wiggler  where  Leff  is  the  effective  interaction 
length.  This  requirement  can  be  expressed  by  the  following  inequality  [2] 

A EWJE  <  4(1  +  4,/2)/(4Z.tf),  (2a) 

where  the  total  longitudinal  energy  spread  is 

A Em  IE  =  s!/(TbE)  +  slL  !(2Rl)  +  n1  K 24)  •  (2b) 

In  Eq.  (2b),  the  first  term  represents  the  fractional  energy  spread  due  to  longitudinal 
emittance  ez ,  the  second  term  arises  from  transverse  rms  emittance  sn  L ,  and  the  third 

term  is  due  to  transverse  wiggler  gradients.  The  right  hand  side  of  Eq.  (2a)  is  the 
approximate  FEL  conversion  efficiency,  where  the  effective  interaction  length,  Lejf  ,  is 

the  gain  length  in  the  case  of  the  amplifier  or  the  wiggler  length  in  the  case  of  the 
oscillator.  If  the  inequality  in  Eq.  (2a)  is  well  satisfied  the  spread  in  the  axial  velocity  of 
the  electrons  plays  little  role  in  the  FEL  interaction.  As  an  example,  for  a  MW-class  FEL 
operating  in  the  IR  regime,  the  typical  parameters  are  snl  ~  15  mm  mrad , 

s:  «  100  keV psec,  Rb  »  0.3mm,  rb  «  lpsec(rms),  E  =  80MeV,  a2w  =  2 ,  and 

=  3  cm  [16J.  For  these  parameters  sz  l(jb  E)<  0.13% ,  e~  L  /(27?2)  »  0.13%, 

iz2  a2wRl  /( 2A2  )  «  0.1%  and  AEtotal/ E  «  0.36% .  For  an  amplifier  having  a  gain  length  of 
30  cm  the  right  hand  side  of  Eq.  (2a)  is  »  5%  and  the  inequality  is  well  satisfied.  For  an 
oscillator  the  right  hand  side  of  Eq.  (2a)  is  -  1/2 Nw ,  where  Nw  is  the  number  of  wiggler 
periods.  For  Nw  =  100  the  inequality  is  also  satisfied. 

d)  Coherent  and  Incoherent  X-Ray  FELs 
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The  FEL  can  generate  coherent,  polarized,  short  pulses  of  x-rays  for  numerous 
applications  in  research.  A  number  of  large-scale  electron  accelerator  facilities 
throughout  the  world  are  being  used  (or  will  be  commissioned  soon)  for  coherent  x-ray 
generation  using  a  conventional  FEL  configuration  [17-20].  Generation  of  x-rays  at  these 
facilities  typically  requires  electron  beam  energies  in  the  multi-GeV  range  with  peak 
currents  in  the  multi-kA  range,  and  wiggler  lengths  of  many  tens  of  meters.  An  x-ray 
FEL  amplifier  can  be  operated  in  the  self-amplified  regime,  eliminating  the  need  for  a 
coherent  input  x-ray  source  [21-23]. 

The  wiggler  field  in  the  FEL  can  be  replaced  with  an  intense  laser  beam 
propagating  anti-parallel  to  the  electron  beam  [24],  In  a  laser  pumped  x-ray  FEL  electron 
beam  energies  in  the  multi-MeV  range  would  be  sufficient,  leading  to  a  more  compact 
device.  However,  a  number  of  major  challenges  would  have  to  be  overcome  before  a 
laser  pumped  coherent  x-ray  FEL  would  be  feasible,  e.g.,  ultra-high  electron  beam 
brightness.  A  compact  source  of  incoherent  x-rays  can  be  configured  in  which 
microwaves  or  laser  radiation  is  Compton  backscattered  off  a  relativistic  electron  beam 
[25-27], 

e)  Backward  Raman  Amplification 

In  the  backward  Raman  amplification  (BRA)  scheme,  two  counter  propagating 
laser  beams  are  resonantly  coupled  by  a  plasma  wave  [28-30].  The  laser  beams  have  a 
frequency  difference  which  is  approximately  equal  to  the  plasma  frequency.  The  higher 
frequency  ( coa ),  long  pulse  pump  laser  supplies  energy  to  the  lower  frequency 

( o\  »  co0  —  cop ),  short  pulse  seed  laser  which  undergoes  amplification  and  pulse 

compression.  The  phase  velocity  of  the  plasma  wave  is  in  the  same  direction  as  the  pump 
laser. 

The  BRA  in  the  linear  regime  is  fundamentally  the  same  as  the  FEL  mechanism 
in  the  frame  of  the  electron  beam.  In  the  electron  beam  frame  the  wiggler  field  is 
essentially  a  counter  propagating  electromagnetic  wave.  Thus,  replacing  the  wiggler 
field  in  Eq.(l)  with  a  counter  propagating  laser  beam  results  in  the  BRA  dispersion 
relation.  By  setting  vz  =  0  in  the  FEL  dispersion  relation,  Eq.(l),  with  appropriate 
redefinition  of  symbols  the  BRA  dispersion  relation  is  obtained 
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( k 2  —co2  / c2 )  ((co  -  co0)2  —  co2)  =  —  2 co2  a20kl  and  sketched  in  Fig.  5.  In  the  BRA 

process  the  seed  laser  grows  and  is  temporally  compressed  at  the  expense  of  the  pump 
laser.  This  process  can  generate  high  intensity  ultra-short  laser  pulses.  Experiments  have 
demonstrated  amplification  of  seed  intensities  by  a  factor  of  104  with  an  efficiency  of  > 
6%  [31,32]. 


CO 


Figure  5.  Dispersion  diagram  describing  the  backward  Raman  amplification  process.  This 
dispersion  diagram  also  describes  the  EEL  process  in  the  beam  frame. 


III.  Ultra  Short  Pulse  Laser  Interactions 

A  great  deal  of  progress  has  been  made  in  the  development  of  ultra-short  pulse 
lasers  (USPL)  [33].  Table  I  lists  some  of  the  important  parameters  associated  with  these 
lasers. 


Peak  power 

1012-  1015W 

Intensity 

1013  -  1023  W/cm2 

Pulse  duration 

10-12  -  10'14  sec 

Average  power 

<  100  W 

Table  I.  Range  of  parameters  for 
ultrashort  pulse  lasers 
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Ultrashort  pulse  lasers  have  numerous  applications,  including  i)  laser  driven  electron  and 
ion  acceleration  (plasma,  Cherenkov,  vacuum);  ii)  laser  synchrotron  (Compton 
backscattering)  source  of  x-rays;  iii)  remote  sensing  and  radioactivity  detection  and  iv) 
remote  atmospheric  lasing. 

a)  Atmospheric  Propagation  of  Ultrashort  Pulse  Lasers 

Atmospheric  propagation  of  USPLs  with  ~  TW  power  levels  and  pulse  durations 
of  ~  100  fsec  is  strongly  affected  by  various,  interrelated  linear  and  nonlinear  processes 
[33-41].  These  include  diffraction,  Kerr  focusing,  group  velocity  dispersion,  spectral 
broadening,  self-phase  modulation,  photo-ionization,  plasma  defocusing,  and  energy 
depletion  due  to  both  ionization  and  ohmic  heating.  A  laser  pulse  propagating  in  air  can 
be  longitudinally  and  transversely  focused  at  remote  distances  (>  km)  to  reach  intensities 
which  will  ionize  the  air.  Longitudinal  pulse  compression  is  achieved  by  introducing  a 
negative  frequency  chirp  on  the  pulse  while  nonlinear  transverse  focusing  is  caused  by 
the  Kerr  effect  (see  Fig.  6).  Asa  result  of  the  focusing,  the  laser  intensity  is  greatly 
enhanced  (>  1013  W/cm2),  resulting  in  tunneling  ionization  of  air  and  a  localized  plasma 
filament  ~  1  m  in  length,  ~  100  pm  in  radius,  having  an  electron  density  of  ~  1016cm“3 
and  a  temperature  of  ~  leV  [40].  While  undergoing  filamentation  (plasma  formation)  the 
intensity  of  the  USPL  is  depleted  due  to  ionization  losses  or  diffraction  and  has  a  typical 
propagation  range  of  many  tens  of  meters.  The  length  of  the  filament  is  determined  by 
air  chemistry  processes  such  as  electron  recombination  and  attachment. 
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Figure  6,  Atmospheric  propagation,  filamentation  and  spectral  broadening 

The  propagation  of  USPLs  in  the  atmosphere  may  have  important  applications  in 
the  areas  of  active  and  passive  remote  sensing,  electronic  countermeasures  and  induced 
electric  discharges,  e.g.,  artificial  lightning.  For  example,  localized  ultraviolet  radiation 
generated  at  a  remote  distance  by  laser  ionization  can  provide  a  source  for  active 
fluorescence  spectroscopy  of  biological  and  chemical  agents  in  the  atmosphere.  Pulses  of 
intense,  directed  white  light  may  also  find  applications  in  the  areas  of  hyperspectral 
imaging  and  differential  absorption  spectroscopy.  The  nonlinearities  associated  with  the 
air  and  plasma  filament  can  result  in  super-continuum  spectral  broadening  of  the  USPL 
and  ultraviolet  radiation. 

Experiments  using  ultrashort,  high  intensity  laser  pulses  have  demonstrated 
atmospheric  propagation,  air  breakdown,  filamentation,  and  white  light  generation  [42- 
45].  Intense,  directed  white  light  pulses  have  been  generated  and  backscattered  from 
atmospheric  aerosols  up  to  altitudes  of  ~15  km  [37].  The  generation  of  THz  radiation  by 
femtosecond  pulses  has  also  been  observed  and  analyzed  [46,47]. 

Filamentation  in  Air 

Hot  spots  on  the  intensity  profile  of  a  laser  beam  can  grow  as  a  result  of  a 
filamentation  instability.  Filamentation,  i.e.,  transverse  break-up  of  a  laser  beam  is  due  to 
the  interplay  between  diffraction  and  nonlinear  self-focusing  (Kerr  effect).  Consider  a 
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laser  beam  propagating  in  a  neutral  gas  for  which  the  nonlinear  focusing  power  is 
PK  =  A2  /(27rnK)  where  nK  is  the  Kerr  effect  contribution  to  the  refractive  index  and  the 
transverse  laser  intensity  profile  is  slightly  perturbed  by  a  small,  localized  hot  spot.  The 
spatial  growth  rate  of  this  filamentation  instability  [40]  is 


T  = 


/l 


n  x. 


3tz1 

v  Pr 


1 


\  1/2 


(3) 


V-L  J 


where  I  is  the  laser  intensity  and  x_|_  is  the  characteristic  transverse  dimension  of  the 
filament,  i.e.,  spot  size.  As  a  function  of  the  dimension  of  the  filament,  the  growth  rate 
vanishes  for  jc±  <  xmin  =  [PK  I(3jt  I)]'12 ,  reaches  a  maximum  equal  to 

rmix  =  3  A I  I(2Pk  )  at  xj_  =V2xmin,  and  decreases  inversely  with  x±  as  x±  increases. 
At  maximum  growth  rate  the  power  within  the  filament  is  roughly  equal  to  PK  .  It  is 
therefore  expected  that  a  laser  beam  with  a  power  P  will  break-up  into  N  filaments  where 
N  <  P  /  PK .  As  an  example,  the  nonlinear  focusing  power  associated  with  air  for  a 
lpm  wavelength,  ~  1  nsec  duration  laser  is  in  the  range  of  ~  3  GW 
( nK  ~5.6xlO-19cm2AV)  [48,49].  For  example,  a  3  TW  laser  pulse  can  eventually  break¬ 


up  into  ~  103  filaments,  see  Fig.7  [42]. 
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Figure  7:  Filamentation  of  a 
400  fs,  A  =  1.05  pmUSPL 
after  13  m  of  propagation  in 
air. 
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The  free  electron  density  in  air  can  change  because  of  ionization,  recombination 
and  attachment  processes.  The  rate  equation  for  electron  density  Ne  is  [50-53] 

^r-  =  WNM-r,N'-p,Nl,  (4) 

ot 

where  W is  the  photoionization  rate,  Nair  =  2.7xl019  cm-3  is  the  air  density,  77  is  the 
electron  attachment  rate  and  (3e  is  the  electron-ion  dissociation  recombination  rate. 
Typical  values  for  the  attachment  and  recombination  coefficients  are  7  -  108  sec-1  and 
/?,  -  3x1 0“8  cm3sec_1 .  For  short  laser  pulses,  free  electrons  are  generated  by  multi¬ 
photon  and  tunneling  processes;  avalanche  ionization  is  not  significant  for  USPLs 

The  length  of  the  plasma  filament  is  determined  by  the  electron  attachment  and 
recombination  coefficients  and  from  Eq.(4)  and  is  given  by 

i  fiiament  ~  c  l(ij  +  j3e  /V,)  <  2m.  The  USPL  loses  energy  through  ionization  and  has 
a  propagation  length  of  Ldepletion  «  P t,  KUwn  Ne  n R]  / 2)  ~  100m,  where  P  is  the 
laser  power,  tl  is  the  pulse  length  and  Uion  is  the  ionization  potential  (see  Fig  8).  The 
balancing  of  Kerr  focusing  and  plasma  defocusing  leads  to  intensity  clamping  of  the 
USPL  which  is  typically  Iclamp  ~  1013-1014  W/cm2. 


filament 

◄ - ► 

- ►  Lprop  •  _  - 

< -  - ► 

Figure  8.  USPL  filamentation,  indicating  filament  length  and  propagation  limits. 

b)  Laser  Driven  Acceleration 

The  extremely  large  field  associated  with  USPLs  propagating  in  plasmas  can 
generate  plasma  waves  which  trap  and  accelerate  electrons  [54-68],  In  the  standard  laser 
wakefield  accelerator  (LWFA),  a  short  laser  pulse,  on  the  order  of  a  plasma  wavelength 
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long,  excites  a  trailing  plasma  wave  that  can  trap  and  accelerate  electrons  to  high  energy 
(see  Fig.9).  Numerous  analyses  and  experiments  have  been  performed  to  characterize 
and  advance  the  LWFA  concept.  However,  there  are  a  number  of  issues  that  must  be 
resolved  before  a  practical  high-energy  accelerator  can  be  developed.  These  include 
Raman,  modulation  and  hose  instabilities  that  can  disrupt  the  acceleration  process  [68]. 

In  addition,  extended  propagation  of  the  laser  pulse  is  necessary  to  achieve  high  electron 
energy. 

In  the  absence  of  optical  guiding  the  acceleration  distance  is  limited  to  a  few 
Rayleigh  ranges  ZR  =  ttR2  /  X  ~  1cm  which  is  far  below  what  is  necessary  to  reach  GeV 

electron  energies  [54].  The  diffraction  limitation  can  be  mitigated  by  employing  a 
plasma  channel  in  which  the  plasma  density  increases  as  r2.  The  physics  of  laser  beams 
propagating  in  plasma  channels  has  been  studied  [56,64,68,69]  and  there  is  ample 
experimental  confirmation  of  extended  guided  propagation  in  plasmas  and  plasma 
channels  [70-76],  In  addition,  dephasing  of  electrons  in  the  wakefield  can  limit  the 
interaction  length  to  (A  /  A)2  Ap  ~  30cm  and,  therefore,  limit  the  energy  gain.  By 

increasing  the  plasma  density  as  a  function  of  distance,  the  phase  velocity  of  the 
accelerating  field  behind  the  laser  pulse  can  be  made  equal  to  the  speed  of  light.  Thus, 
electron  dephasing  in  the  accelerating  wakefield  can  be  postponed  and  energy  gain 
increased  by  spatially  tapering  the  plasma  channel.  Staging  LWFAs  can  be  used  to 
overcome  the  laser  energy  loss  limitation.  Injection,  trapping  and  acceleration  in  the 
wakefield  continues  to  be  of  interest  [77]. 
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Figure  9.  Schematic  of  laser  wakefield  accelerator  mechanism 


A  simplified,  one-dimensional  model  of  the  laser  wakefield  accelerator  driven  by 
the  ponderomotive  force  induced  by  the  USPL  is  shown  in  Fig.9  [54],  The 
ponderomotive  force  Fpond  =  <?(vxB/c)z  =  -  (c/2 /2maj2)d(EL-EL)/d z  associated  with  the 


laser  pulse  generates  large  amplitude  plasma  waves  (  >  100  GV/m)  that  can  accelerate 
self-trapped  or  injected  electrons  to  high  energies. 

The  plasma  density  wave  (wakefield)  is  driven  by  the  ponderomotive  force  and 
satisfies  a  driven  harmonic  oscillator  equation 


d ?SNp  ,  2  d 

- ^  +  ca\SN0  =  N  c2 — 7 

dt2  p  p  po  dz2 


qE 

2mcco 


(5) 


where  co  is  the  plasma  frequency  N  is  the  unperturbed  plasma  density,  5  N  is  the 


perturbed  density  (wakefield)  and  EL(z„t)  =  E(z,t)sm(k  z  —  cot)  the  laser  field. 

The  accelerating  wakefield  (space  charge  field)  is  given  by  dEaccel  I  dz  =  4  nqSN  p. 


The  accelerating  gradient  \Eaccel\  is  proportional  to  ~  p  1 ,  where  I  =  c 


/8 k  is 


the  intensity,  and  can  be  ~  Kr  times  greater  than  accelerating  gradients  in  RF  accelerators. 
i)  Optical  Guiding  in  LWFA 
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The  laser  spot  size  R(z )  as  a  function  of  axial  position  in  a  plasma  with  refractive 


index  n  is  given  by  [64] 


a2  R  _  /L2  2 J  dn\ 

dz2  ~  tt2R 3  +  R\  dr  j 


(6) 


where  (...)  is  an  average  performed  using  the  laser  intensity  as  the  weighting  function. 

The  full  refractive  index  associated  with  the  laser  plasma  interaction  having  a  parabolic 
density  variation  in  the  radial  direction  is 


The  expression  for  the  refractive  index  is  obtained  by  noting  that 
ojp(r)  =  ( 47rq 2  N p(r) / y(r)m)1'2  is  the  plasma  frequency,  where 

Np(r )  =  Npo  +  An  r2  /  R2  is  the  plasma  density  in  the  channel,  y(r)  is  the  relativistic 


mass  factor,  j3os(r)  =  q  E(r) Kmc  co) ,  copo  =a>po{r  =  0)  and  Npo  =Np(r  =  0).  For  short 
laser  pulses  ctl  <  Ap  the  relativistic  term  J3PS  /  4  just  cancels  the  wakefield  term 


SN p  /  Npo .  In  the  short  pulse  limit  the  laser  spot  size  dynamics  in  a  plasma  channel  is 
given  by 


8-J?  _  RU  _ 

5?  ”  l  )-&’ 


where  nc  =  \  ti  n  r,  R~ )  is  the  depth  of  the  plasma  channel  and  r  is  the  classical  electron 

radius  (see  Fig.  10).  Note  that  in  the  absence  of  a  plasma  wave  and  plasma  density 
channel,  i.e.,  uniform  plasma  density,  optical  guiding  is  possible  through  a  relativistic 
effect.  Here  the  laser  power  must  equal  the  relativistic  focusing  power,  which  is 
Prel[GW]  =  l7.4(co/cDp)2. 
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Figure  10.  Guided  propagation  of  Gaussian  laser  beam  in  a  parabolic-profile  plasma 
density.  In  a  plasma  channel  the  guiding  condition  is  independent  of  wavelength. 

Extended  optical  guiding  is  important  for  achieving  significant  energy  gain  in  the  LWFA. 
Laser  wakefield  experiments  at  LBNL  employing  optical  guiding  of  a  40  TW  laser  pulse 
in  a  gas  jet  configuration  resulted  in  observation  ~  GeV  electrons  [78], 

ii)  Raman  and  Modulational  Instabilities 

The  dynamics  of  a  laser  pulse  propagating  in  a  uniform  plasma  channel  is  affected 
by  Raman  and  modulational  instabilities.  In  the  broad  beam  limit  ( Ra  »crL),  however, 

the  nonlinear  terms  in  the  wave  equation  that  lead  to  Raman  and  modulation  instabilities 
tend  to  cancel  [79].  This  can  result  in  pulse  propagation  over  extended  distances,  limited 
only  by  dispersion  with  minimal  pulse  distortion  due  to  instabilities. 

iii)  LWFA  in  a  Tapered  Plasma  Channel 

By  employing  a  tapered  plasma  channel  electron-wakefield  dephasing,  i.e., 
slippage,  can  be  postponed  [80,81],  The  phase  of  the  wakefield  behind  the  laser  pulse  is 
y/(z,z)  =  cop{z)(j  —  rL! 2)  and  the  phase  velocity  is  given  by 


(9) 


where  z  =  t—  I  dz'/vg(z')  is  the  time  measured  in  the  frame  of  the  pulse  with  group 


0 


velocity  \g  (z)  and  zL  is  the  laser  pulse  length.  The  phase  velocity  of  the  wakefield 
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varies  with  distance  behind  the  pulse.  The  phase  velocity  increases  (decreases)  with 
distance  from  behind  the  pulse  for  an  increasing  (decreasing)  plasma  density.  The 
location  behind  the  pulse,  tc  ,  for  which  the  phase  velocity  equals  the  speed  of  light  in 


vacuum  is  given  by 


M  %  /(2&r)  +  2c  l(co20  R20)  Ctl 

CTr(z)  =  — -  +  - ~ 

d\ncop/dz  2 


(10) 


The  solid  curve  in  Fig.  1 1  shows  the  energy  of  a  test  electron  along  a  tapered  plasma 
channel  (increasing  plasma  density)  that  is  trapped  and  accelerated  to  ~  4  GeV.  The 
dashed  curve  shows  an  energy  gain  of  ~  lGeV  in  an  untapered  channel. 


Figure  1 1 .  Energy  gain  as  a  function 
of  propagation  distance  normalized 
to  the  Rayleigh  range  ZR,.  Dashed 
curve  is  for  an  untapered  channel 
(maximum  energy  ~  1  GeV)  while 
the  solid  curve  is  for  a  tapered 
channel  (maximum  energy  ~  4  GeV). 


Spatially  periodic  plasma  channels,  with  periods  as  small  as  ~  50  jum,  have  been 
demonstrated  [82],  The  modulated  plasma  channel  can  support  lower  phase  velocity 
wakefields  which  can  more  readily  trap  low  energy  electrons  and  accelerate  them  into 
higher  phase  velocity  wakefields. 

IV.  Laser  Propagation  in  Atmospheric  Turbulence 

Atmospheric  turbulence,  absorption,  scatterings,  scintillation  and  thermal 
blooming  can  have  an  important  effect  on  laser  propagation  [83].  The  following 
discussion  deals  primarily  with  laser  propagation  through  atmospheric  turbulence.  The 

refractive  index  of  turbulent  air  is  n  =  1  +  dn(x,  y,  z,  t ) ,  where  {Sri)  =  0  and  (S n2')  ^  0 . 
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The  strength  of  the  turbulence  is  characterized  by  the  structure  parameter  C2  which 
typically  is  in  the  range  of  1CT15  —  1  (T13  m-2/3  at  ground  level  [83-85]. 

Turbulence  leads  to  an  increase  in  spreading  of  the  laser  beam  spot  size,  Rs .  In 
addition,  turbulence  leads  to  wandering  of  the  laser  beam  centroid,  Rw  Figure  12  shows 

the  laser  beam  spot  (small  circles)  at  several  instants  in  time.  Over  a  time  scale  that  is 
long  compared  with  that  associated  with  wander  one  obtains  the  broadened  spot  shown 
by  the  large  dotted  circle  in  Fig.  12. 

The  time  averaged  laser  spot  size  R(L )  at  range  L  is  given  by 

R\L )  =  (A,  LI  n  R0)2  (M4  +  2.9  {RJr0)2)  +  G2jitterL2  +  R2(  1  -  L/Lf)2,  (11) 
where  R~(R;v  +  R2)m ,  Ra  is  the  initial  spot  size,  A  is  the  laser  wavelength,  M 2  is  the 
intrinsic  laser  beam  quality,  Lf  is  the  focal  length,  r0  is  the  transverse  coherence  length 
associated  with  turbulence  and  0jitter  is  the  mechanical  jitter  angle  of  the  beam.  For  a 

constant  value  of  C2  the  transverse  coherence  length  is  ra  =  0.184(/l2  /C2  L)3/5 . 
However,  as  far  as  laser  beam  propagation  is  concerned  the  effective  strength  of 
turbulence  is  not  determined  by  C2  alone,  but  also  depends  on  the  range  of  propagation 
and  the  wavelength.  The  relevant  combination  is  the  Rytov  variance,  defined  by 
<j\  =  10.5  C2  A~116  L11'6 .  For  weak  turbulence  ( a\  <  0.3),  beam  wander  and  spreading 
are  distinct  and  tip-tilt  compensation  can  be  used  to  correct  for  beam  centroid  wander. 

On  the  other  hand,  for  strong  turbulence  ( o\  >  0.3)  the  beam  breaks  up  into  multiple 
beams  and  beam  wander  and  spreading  are  not  distinct,  hence,  tip-tilt  compensation  is 
less  effective.  The  Rytov  variance  is  a  measure  of  the  intensity  scintillation  level  on  axis 

at  the  target,  i.e.,  g\  =  [{l  —  (/))2  ^/(/)2  for  weak  turbulence. 

Figure  12.  Schematic  showing  the  effects  of 
turbulence  on  beam  propagation,  Rw  is  the 

centroid  displacement  (wander),  and  Rs  is  the 
increase  in  spot  size  (spreading),  where 
R*(R2 +R2)1/2. 

X 
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Figure  13  shows  the  laser  intensity  contours  for  three  levels  of  atmospheric  turbulence. 
In  limit  of  extremely  weak  turbulence  ( C2  =  0  ,  ct2r  =  0 ,  Fig.  13(a))  the  laser  intensity  is 
well-defined  with  a  relatively  small  spot  size  [86],  In  moderate  turbulence 
(C2n  =  10 ~14 m~2/3  ,<yR  »  0.1,  Fig.l3(b))  the  laser  profile  is  distorted  and  the  spot  size  is 

larger.  In  strong  turbulence  ( C2n  =  10-13  m~2n  ,  <j2r  »  1 ,  Fig.  13(c))  the  laser  profile  is 

highly  distorted.  For  sufficiently  strong  turbulence  a  single  discernable  high-intensity 
region  in  the  cross-section  may  not  be  obtained. 
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Figure  13.  Comparison  of  the  intensity  contours  in  the  x-y  plane  at  a  range  of  0.5  km  for 
three  levels  of  turbulence  using  the  propagation  code  HELCAP  [86]. 
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Effective  beam  control  of  high  energy  lasers  (HELs)  is  a  key  component  of  a  laser 
directed  energy  system.  Strong  turbulence  can  have  a  significant  deleterious  effect  on  the 
propagation  of  the  HEL  beam.  By  appropriately  modifying  the  wavefront  of  the 
transmitted  beam  the  laser  power  can  be  more  effectively  focused  on  target  [87,88],  A 
beacon  beam  can  be  used  to  record  the  atmospheric  phase  aberrations  used  to  modify  the 
amplitude  and  phase  of  the  transmitted  high  power  beam;  see  Fig.  14  for  a  simplified  view 
of  the  process.  The  information  needed  for  modifying  the  amplitude  and  phase  is 
obtained  by  phase  conjugating  the  beacon  beam,  i.e.,  exp(i k  •  r)  =>  exp(— i k  •  r) .  Phase 
conjugation  can  be  achieved  by  using  deformable  mirrors  or  by  a  nonlinear  optical 
mixing  mechanism,  e.g.,  4-wave  Brillouin  mixing  [89,90]. 
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Figure  14.  Using  phase  and  amplitude  information  provided  by  a  beacon  beam  the 
distortion  of  the  transmitted  high  power  laser  beam  is  gradually  undone  as  it  propagates 
to  the  target. 

V.  Applications  of  High-Power,  High-Intensity  Lasers 

Advances  in  FELs  and  solid-state  lasers  (including  fiber  lasers)  have  made  them 
prime  candidates  for  high-average  power  applications.  These  include  power  beaming  and 
directed  energy  applications.  In  this  section  we  discuss  remote  detection  of  radioactive 
material  and  remote  atmospheric  lasing  for  chemical/biological  detection,  employing 
high-intensity  lasers. 

a)  Remote  Detection  of  Radioactive  Material 
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A  recently  proposed  radioactivity  detection  concept  is  based  on  a  high  power  THz 
pulse  inducing  avalanche  breakdown  and  spark  formation  in  the  vicinity  of  the 
radioactive  material  [91,92],  In  this  detection  concept  the  time  delay  for  spark  formation 
and  the  breakdown  rate,  which  are  functions  of  the  initial  ion  and  electron  densities,  can 
provide  a  direct  signature  for  the  presence  of  radioactivity. 

We  discuss  a  new  concept  which  is  based  on  the  detection  of  electromagnetic 
signatures  in  the  vicinity  of  radioactive  material  and  can  enable  stand-off  detection  [93], 
Radioactive  materials  emit  gamma  rays  which  ionize  the  surrounding  air.  The  ionized 

electrons  rapidly  attach  to  oxygen  molecules  forming  0~  ions.  The  density  of  0~ 
around  radioactive  material  can  be  several  orders  of  magnitude  greater  than  background 
levels.  The  elevated  population  of  0~2  extends  several  meters  around  the  radioactive 
material.  Electrons  are  easily  photo-detached  from  0~  ions  by  laser  radiation.  The 
photo-detached  electrons,  in  the  presence  of  laser  radiation,  initiate  avalanche  ionization 
which  results  in  a  rapid  increase  in  electron  density.  The  rise  in  electron  density  induces 
a  frequency  modulation  on  a  probe  beam  which  becomes  a  direct  spectral  signature  for 
the  presence  of  radioactive  material  [93], 

High-power,  short-pulse  lasers  propagating  in  the  atmosphere  can  come  to  a  focus 
in  the  vicinity  of  the  radioactive  material  by  making  use  of  longitudinal  compression  and 
transverse  focusing,  see  Sec.  Ilia  and  [40].  The  negative  ions  produced  by  the 
radioactive  material  have  an  electron  affinity  resulting  in  a  low  ionization  potential  of 
0.46  eV  and  can  be  photo-detached  by  laser  radiation  ( ~  0.8-lpm ).  The  detection 
concept  is  based  on  a  probe  radiation  beam  undergoing  a  frequency  modulation  while 
propagating  in  a  temporally  increasing  electron  density.  The  frequency  modulation  on 
the  probe  beam  becomes  a  spectral  signature  for  the  presence  of  radioactive  material.  A 
schematic  of  the  detection  concept  is  shown  in  Fig.  15. 


22 


Photo-detaching 
and  ionizing 
laser  beam,  <d 


Enhanced 
levels 
of  02- 


Shielding 

Radioactive 

material 


Incident  Probe  Frequency ,  coB 
Reflected  Frequency,  copmbe(t )  =  a>0  +  Aoj(t) 


Elevated 

electron 

density 


Ionizing 

gamma 

radiation 


Figure  15.  Schematic  of  active  remote  radioactivity  detection  concept.  Laser  radiation 
photo-detaches  electrons  from  O”  ions,  providing  electrons  for  an  avalanche  (collisional) 
ionization  process  which  increases  the  electron  density.  A  probe  beam  in  the  presence  of 
a  temporally  increasing  electron  density  undergoes  a  frequency  modulation  which 
becomes  a  spectral  signature  of  radioactivity. 


To  obtain  the  frequency  modulation  on  a  probe  pulse  it  is  necessary  to  follow  the 
time  evolution  of  the  electron  and  negative  ion  density,  which  are  sensitive  functions  of 
air  chemistry  processes  and  electron  heating  by  the  laser  radiation.  The  source  terms  for 
the  electrons  include  radioactivity,  detachment,  photo-detachment  and  photo-ionization, 
while  the  loss  terms  include  various  attachment  and  recombination  processes  including 
aerosols.  The  expressions  for  the  rate  of  change  of  electron  density  N e  and  negative  ion 


density  N_,  i.e.,  02  in  the  present  context,  are  [93], 

dNJdt  =  (l+amd)Qrad  +  Se-  Le,  (12a) 

dN_/dt  =  S_  —  L_,  (12b) 

where  Se  represents  the  various  electron  source  terms,  Le  is  the  electron  loss  terms,  S_ 
represents  the  ion  sources,  L_  is  the  ion  loss  terms  and 

VradKr  ( Ee )  1  eXp  R I  Ly) 


a 


rad 


47r(Lr_e)  A E  Q] 


' rad 


R 


(13) 


is  the  radiation  enhancement  factor.  For  a  detailed  discussion  of  the  various  terms  in 
Eq.(12a,b)  see  Ref  [93], 
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Since  the  free  electron  density  is  generally  much  less  than  the  ion  density, 

Ne  «  N_  the  steady  state  electron  and  negative  ion  densities  are  given  by 

Ne  *  iPn  Nn '  7> +  ■ aj)  OZ  /  J+  ,  and  N_  ~^(i  +  anill)Qrad//3+  ,  where 
J3+  ~  2xl0-6  cm3  /  sec  is  the  recombination  rate,  rj  «  108  sec-1  is  the  attachment  rate 
and  J3n  ~  (5  - 10)  x  10  19  cm3  /sec  is  the  negative  ion  detachment  rate  due  to  collisions 
with  neutrals.  In  the  absence  of  radioactive  material  ( amd  =  0 )  the  ambient  (background) 
electron  and  negative  ion  density  are  Ne  ~  10-3  cm-3  and  N_  «  3xl03  cm-3 ,  where 
Qrad  =  20  cm-3  sec-1 .  At  a  distance  of  4  m  from  a  radioactive  source  containing  10  mg 
of  60 Co ,  the  radioactive  enhancement  factor  is  arad&2xl04  and  the  elevated  electron 

and  negative  ion  densities  are  Ne  »  0.2cm-3  and  N_  «  4.7xl05  cm-3 . 

The  various  electron  and  ion  source  and  loss  terms  in  Eqs.(12a)  and  (12b),  in 
particular  the  collisional  ionization  rate,  are  functions  of  the  electron  temperature.  The 
electron  temperature  is  determined  by  the  collisional  electron  heating  (Ohmic  heating)  by 
the  laser  radiation  and  the  cooling  effect  resulting  from  excitation  of  vibrational  modes  of 
the  air  molecules.  The  equations  for  the  electron  temperature  Te  is  [93], 

(3/2)  3(/V,  r,)/0(=  (J-E)  -  (H2)(NelTcocl)(Tt -TeJ,  (14) 

where  Te  is  the  electron  temperature,  (j  -E)  is  the  Ohmic  heating  rate,  rcool  is  the 
electron  cooling  time  due  to  inelastic  collisions  and  Teo  =  0.025  eV  is  the  ambient 
electron  (room)  temperature. 

A  probe  beam  with  frequency  co0  propagating  through  a  region  of  space  in 

which  the  electron  density  is  changing  with  time  will  undergo  a  frequency  change.  The 
electron  density  in  the  vicinity  of  the  radioactive  source  and  under  the  influence  of  the 
laser  radiation  varies  in  space  and  in  time.  The  frequency  of  an  electromagnetic  probe 
beam  propagating  in  such  plasma  will  vary  in  space  and  in  time.  As  an  illustration  we 
consider  the  case  where  the  rise  in  electron  density  is  exponential  in  time  and  spatially 
uniform  within  a  region  0  <  z  <T.  Here  the  plasma  frequency  is  given  by 

co2p(z,t )  =  O)2poexp(viont) g(z) ,  where  vion  is  the  effective  ionization  rate  and 
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g(z)=l,  0  <  z  <L,  defines  the  axial  extent  of  the  plasma  region.  The  frequency 
modulation  on  the  probe  beam  is  given  by  coprobe(z,t)  =  co0  +  A co(z,t),  where 

A<x>(z,t)  =  co2po  1(2 coa) exp (vion 0(1  -  exp(-i/.„,  z/c)).  (15) 

The  maximum  fractional  frequency  shift  occurs  for  z  >  c  /  vion  ~  1  cm  and  is 
Aojma  /  coa  =  (a>po  /  2&r )  exp  ( vion  t ) .  The  effective  ionization  rate  can  vary  widely  but  is 

typically  vion  ~  1011  sec-1 .  The  frequency  modulation  on  the  probe  beam  can  be 

substantial  and  is  readily  measurable. 

To  illustrate  the  radioactivity  detection  concept  we  consider  the  following 
example.  In  the  numerical  examples  the  system  quantities,  i.e.,  electron  density,  ion 
density,  electron  temperature  and  frequency  modulation  on  a  probe  laser,  are  calculated 
both  in  the  absence  ( arad  =  0 )  and  in  the  presence  of  radioactive  material  by  numerically 
solving  the  full  set  of  coupled  equations.  The  radioactive  source  is  assumed  strong 
enough  to  produce  a  radiation  enhancement  factor  of  arad  =  103 .  We  take  the  ionizing 

laser  to  have  a  peak  intensity  of  I peak  =  160  GW/cm2  and  pulse  duration  of  zlaser  =  lnsec . 
In  these  examples,  the  probe  beam  is  taken  to  be  a  millimeter  wave  source  of  frequency 
f  probe =  94  GHz ,  ( Aprobe  =  3.2mm).  The  background  radiation  is  taken  to  be 

Qmd  =  30  disintegrations  /(cm 3  -sec) . 
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Figure  16.  Electron  density  as  a  function  of  time,  (a)  background  radiation  only 
( amd  =0),  (b)  in  the  presence  of  radioactivity  ( amd  =  103).  The  laser  parameters  are, 

/l  =  lpm,  I peak  =  160  GW/cm2  and  zlaser  =  lnsec. 

In  the  absence  of  a  radioactive  source,  Fig.  16(a),  the  ionizing  laser  intensity  is  just 
below  the  breakdown  level,  i.e.,  the  electron  density  is  low,  and  there  is  virtually  no 
frequency  modulation  on  the  probe  beam.  Figure  16(b)  shows  the  electron  density  as  a 
function  of  time  in  the  presence  of  radioactive  material  ( arad  =  103 ).  The  electron 
density  at  the  end  of  the  ionizing  laser  pulse  approaches  the  value  of 
Ne  =  1013  cm-3  which  is  an  order  of  magnitude  less  than  the  critical  electron  density. 

The  frequency  modulation  on  the  probe  millimeter  wave  beam  is  shown  in  Fig.  17. 
In  the  absence  of  radioactive  material  there  is  no  frequency  modulation  on  the  probe. 
However,  for  amd  =  103  the  fractional  frequency  modulation  is  significant  and  equal 
to  -  5%,  which  is  readily  detectable. 


Figure  17.  Fractional  frequency  shift 
Ao)/ coa  [%]  versus  time  in  the  presence  of 

radioactive  material  arad  =  103  at  the  probe 
interaction  distance  of  L  =  10cm.  The  laser 
parameters  are  the  same  as  in  Fig.  16. 
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b)  Atmospheric  Lasing 

Nitrogen  lasers  are  typically  based  on  collisional  excitation  using 
electrical  discharges  [94-96].  In  this  section  we  discuss  a  remote  atmospheric  lasing 
mechanism  [97]  in  which  an  ultrashort  pulse  laser  forms  a  plasma  filament  (seed 
electrons)  by  tunneling  ionization  and  a  heater  pulse  thermalizes  the  seed  electrons.  The 
thermal  electrons  collisionally  excite  nitrogen  molecules  and  induce  lasing  in  the 
ultraviolet.  A  remote  atmospheric  lasing  configuration  based  on  collisionally  exciting  the 
N2  lasing  line  at  X  =  337  nm  in  a  heated  plasma  filament,  as  illustrated  in  Fig.  18.  The 
lasing  gain  is  sufficiently  high  to  reach  saturation  within  the  length  of  the  plasma  filament. 
A  remotely  generated  ultraviolet  source  may  have  applications  for  standoff  detection  of 
biological  and  chemical  agents. 


Heater  pulse  maintains  and 

thermalizes  seed  electrons  UV  lasing  u£>pL 


filament  (seed  electrons) 
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2 
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Figure  18.  Schematic  diagram  of  remote  lasing  configuration.  An  ultrashort  pulse 
laser  creates  a  plasma  filament  of  seed  electrons  that  is  heated  by  a  secondary 
heater  pulse.  In  the  energy  level  diagram  for  N2  (right),  the  energetic  electrons 
collisionally  excite  the  N2  molecules  and  induce  lasing. 


A  heater  laser  pulse  is  used  to  maintain  and  thermalize  the  electrons.  Thermal 
electrons  of  a  few  eV  can  induce  a  population  inversion  in  N2  leading  to  the  vibrational- 

electronic  transition  C3Ilu  — >  B3T I  at  337  nm.  In  the  present  configuration  the  UV 

radiation  is  generated  in  both  the  forward  and  backward  directions  along  a  filament  with 
a  gain  sufficiently  high  to  reach  saturation.  For  a  1pm  wavelength  heater  beam  the  Kerr 
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nonlinearity  can  be  the  source  of  third  harmonic  generation  and  a  seed  for  lasing  at  337 
nm. 

To  analyze  the  lasing  process  the  density  matrix  and  Maxwell  equations  [97] 
together  with  the  electron  heating  and  ionization  equations  are  solved.  The  N2  lasing 
model  consists  of  an  open  two  level  system  denoted  by  levels  3  (upper)  and  2  (lower) 
respectively.  Since  levels  3  and  2  are  weakly  excited  the  population  of  level  1  (ground) 
is  taken  to  be  fixed.  The  population  of  the  excited  levels  are  given  by 

d  Ni  I  d  t  =  vCE  n  N  +  vCe,i 3^2  —  Fs  ^3  —  ^stim^N ,  (16a) 

d  N 2 1  d  t  =  vCE  j  2  N  +  (yCD  32  +  r32)^3  —  v2  N  2  T  vstim&N ,  (16b) 

where  r  =  t  —  z  /  c ,  z  is  the  position  within  the  filament,  Nn ,  are  densities  of  the  nth 
level,  TV  »  iVj  =  0.8 Nair  is  the  ground  state  nitrogen  density,  Nair  =  2.7xl019cm-3  is  the 
air  density,  v3  =  vCD  31  +  vCD32  +  T31  +  T32  is  the  decay  rate  of  level  3, 

v2  =  v cd. 2 1  +  vce,23  +^21  is  the  decay  rate  of  level  2,  r.  ;  is  the  radiative  (spontaneous) 

decay  from  level  i  to  j,  vCE  lJ  is  the  collisional  excitation  rate  from  level  i  to  j,  and  vCD  iJ 
is  the  collisional  de-excitation  rate  from  level  i  to  j  (i,  j  =1,2,  3),  vstim  =  ostim  I  Itico  is  the 
stimulated  emission  rate,  I  is  the  UV  laser  intensity,  co  is  the  lasing  frequency  <Jstim  is 
the  stimulated  emission  cross  section  and  AN  =  N3  —  N2.  The  collisional  excitation 

rates  are  dominated  by  electron  excitations  while  the  de-excitation  rates  are  dominated  by 
molecular  collisions. 

The  lasing  intensity  is  given  by 

(8/Sz  +  rj/  =  «/„„/(«■„,)  +  <rjnm  AN  1 ,  (17) 

where  Td  is  the  spatial  damping  rate,  e  is  a  geometric  filling  factor  associated  with  the 
seed  radiation  within  the  plasma  filament,  I stim  =hcocN3  and  zrad  is  the  radiative 

lifetime  from  level  3  to  2  including  the  effects  of  collisions. 

Figure  19  plots  the  temporal  profile  of  the  UV  radiation  and  population  inversion 
at  various  positions  along  the  filament.  Here,  the  heating  pulse  duration  is  ra  =  0.4  nsec , 

and  the  intensity  is  I0  =  7x10' 1  W/cm2 .  Figure  19  shows  a  UV  pulse  (~  0.1  nsec  long) 
forming  near  the  back  of  the  inversion  region.  As  the  UV  pulse  grows  AN  is  depleted 
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and  the  peak  of  the  UV  pulse  moves  forward  relative  to  the  heater  pulse.  After  ~  2  cm  of 
propagation,  the  peak  UV  intensity  is  ~  8  MW/cm2  and  the  fluence  is  ~  0.5  mJ/cm2.  The 
peak  electron  temperature  is  typically  ~  2.5  eV  within  the  heater  pulse. 

Lasing  will  occur  in  both  the  forward  and  backward  directions  when  the  length  of 
the  heated  plasma  filament  is  longer  than  the  saturation  distance.  For  the  numerical 
examples  given  here  this  condition  is  satisfied.  The  length  of  the  heated  plasma  filament 
is  ~  6  cm,  while  the  lasing  saturation  lengths  are  typically  ~  1-3  cm.  It  should  be 
mentioned  that  other  molecular  transitions  may  also  be  excited  and  result  in  lasing,  e.g., 
N2  and  02 . 
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Figure  19.  UV  intensity,  I  (solid 
curves)  and  normalized  population 
inversion  A N / N  =  (N3-N2)/ N 
(dashed  curves)  as  functions  of  time 
at  (a)  z  =  0.9  cm,  (b)  z  =  1.4  cm, 
and  (c)  z  =  2  cm,  for  a  heating  pulse 
with  =  ljum,  I0z0  =  285J/cm2 
and  pulse  duration  x0  =  0.4  nsec 


VI.  Concluding  Remarks:  A  number  of  topics  associated  with  high-power,  high- 
intensity  lasers  have  been  discussed.  The  aim  has  been  to  briefly  illustrate  the  varied 
physical  processes  and  phenomena  that  are  manifested  in  the  interaction  of  these  lasers  with 
matter.  We  have  also  described  some  of  the  actual  and  potential  applications  of  high-power 
as  well  as  high-intensity  lasers.  This  paper  is  not  intended  to  be  an  exhaustive  overview  of 
the  various  topics;  important  references  on  the  topics  may  have  been  omitted. 
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